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Abstract 
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it was shown that at four loops, the anomalous dimensions derived 
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this paper we obtain the remaining wrapping part of these anomalous 
dimensions from string theory and show that these contributions ex- 
actly cure the problem and lead to agreement with both LO and NLO 
BFKL expectations. 



*e-mail: bajnok@elte.hu 

^ e-mail: uf r j anikOif . uj . edu . pi 

^e-mail: tomaszlukowski@gmail . com 



1 



1 Introduction 



The integrable structures which appear both on the gauge theory side [H [2j 
[3j HI [5] and the string theory side [6] of the AdS / CFT correspondence [7] 
give hope to find the full spectrum of both theories. On the gauge theory 
side we are interested in finding the anomalous dimensions of gauge invariant 
operators in M = 4 super Yang-Mills theory, while on the string theory side 
we would like to determine the energies of the superstring excitations in the 
AdS$ x S 5 background. 

A lot of progress has been done for specific type of operators, namely with 
large numbers of fields, and for the corresponding infinitely long/ fast strings 
(strings with large charges). The S-matrix for elementary excitations has 
been predicted from the symmetry algebra up to the overall scalar factor [8] 
which was finally fixed in [9l [10] using the crossing symmetry [TT] . The Bethe 
Ansatz Equations have been derived ^[TSlH^mHSllISlinillBl § ivin g the 
spectrum of states with large quantum numbers. 

Nevertheless, it is known that the Bethe Ansatz is not valid for short 
operators and strings with small charges. On the gauge theory side it is 
due to the fact that at the order of X L , where L is the length of the gauge 
invariant operator, wrapping corrections start to play a role. In an analo- 
gous manner [T3], on the string theory side virtual corrections around the 
worldsheet cylinder appear. In [20] the leading exponential finite size effects 
at strong coupling for a single giant magnon were computed from the iden- 
tification with Liischer-like corrections (see also [211 1221 [23]) and turned out 
to agree exactly with the expression obtained from classical string solution 
[21]. At weak coupling, wrapping corrections were analyzed from the gauge 
theory perturbative point of view in [25]- The wrapping contribution to 
the Konishi operator anomalous dimension at four loops has been explicitly 
computed within perturbative gauge theory [26] (a subsequent independent 
ab-initio perturbative computation of [27J which included both the wrapping 
and non- wrapping graphs reaffirmed this result). The same result was ob- 
tained purely within string theory from the AdS$ x S 5 string sigma model in 

m. 

The wrapping effects appear in the most manifest way for the shortest 
possible operator. The simplest non protected operators, which are general- 
izations of the Konishi operator, are twist two operators: tr ZD 2M Z. Apart 
from being a natural testing ground for wrapping effects, which should ap- 
pear at 4 loops, their anomalous dimensions are also interesting for their 
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own sake as they are intertwined with various seemingly unrelated physical 
processes and have rich analytical properties. Moreover these operators (in 
fact their very close analogues) are also of big importance in ordinary QCD 
(see in particular [29] in the context of the present paper). 

Their anomalous dimensions have been computed perturbatively (in QCD) 
up to 3 loops [30] and using the maximum transcendentality conjecture [21] 
for M = 4 SYM the corresponding answer in M = 4 SYM has been ex- 
tracted [32j [33j [31] . This has been found to agree with the prediction from 
the asymptotic Bethe Ansatz. Subsequently the Bethe Ansatz answer for 
four loops was computed |35j and found to satisfy the maximal transcen- 
dentality principle. It was, however, demonstrated that the Bethe Ansatz 
answer is in conflict with predictions from the BFKL equation [36J. It was 
pointed out that the missing part of the four-loop answer should come from 
wrapping effects. The motivation of our paper is to fill this gap and find the 
leading wrapping correction for general twist two operators and reexamine 
the consistency with BFKL. 

Thus our aim is to extend the calculations of the wrapping corrections to 
the case of arbitrary twist two operators. We will adapt the method used in 
the case of Konishi operator in [28], based on the integrability properties of 
the worldsheet quantum field theory of the superstrings in AdS$ x S 5 . The 
leading contribution from wrapping will appear at order A 4 which exactly 
matches the predictions coming from TBA type considerations [19J. 

The operators considered in this paper are part of a larger family of 
interesting operators from the si (2) sector made up from M derivatives 
(spin) and J complex scalars (twist): tr (D Sl Z . . . D SJ Z). For high er twist, 
there is also a wide range of interesting phenomena e.g. the large spin (M) 
and twist (J) behaviour of their anomalous dimension A(M, J) is determined 
by the BES/FRS integral equations [TUl I3"T] . Remarkably, in the limit when 
both M, J — > oo such that j = is kept fixed the leading logarithmic 

behaviour is governed by the cusp anomalous dimension and the ground-state 
energy density of the 0(6) a model [38l [391 SOI SH S21 S3J SH SSI HSl SZ] . 
We expect that the methods presented in the present paper should also be 
applicable for the higher twist case. 

The plan of this paper is as follows. In section 2 we briefly introduce the 
twist two operators and their main properties. In section 3 the formula for 
wrapping correction for general twist two operator is given. In section 4 we 
point out all ingredients needed to compute the leading wrapping correction. 
All relevant calculations ispresented in section 5, while the final answer is 
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obtained in section 6. Section 7 shows that the wrapping corrections do 
not change the large M asymptotic behaviour of anomalous dimension for 
twist two operators, consequently, the cusp anomalous dimension ('scaling 
function') is untouched by wrapping. In section 8 we analytically continue the 
final result to M = — 1 + uj and then compare it with the BFKL prediction. 
We close the paper with a discussion and several appendices. 

2 Twist two operators and BFKL 

Twist two operators in M = 4 SYM are operators which are made from two 
Z fields (J = 2) and an arbitrary (but even) number of derivatives M in 
a fixed light cone direction. Schematically these operators are thus of the 
form tr ZD M Z. For each even M there is a unique highest weight non BPS 
operator and so we may define the anomalous dimension 

oo 

A(M, J = 2) = 2 + M + ^ MM)g 21 (1) 

i=i 

where g 2 = X/16tt 2 . Here the index I denotes the perturbative loop order. 
This quantity can be computed from the asymptotic Bethe ansatz in the 
si (2) sector exactly up to 3 loops. The answer from the Bethe ansatz at 4 
loops and higher will have to be supplanted by the contribution of so-called 
'wrapping interactions'. The aim of this paper is to compute this contribu- 
tion at 4 loop order from the string sigma model in AdS§ x S 5 . This arises 
due to the identification of the anomalous dimensions with energies of string 
states in AdS$ x S 5 . Since the worldsheet QFT is integrable, and we know 
the exact S-matrix, we know the full on-shell data of the worldsheet QFT 
at infinite volume. We may therefore study the spectrum of energies around 
the infinite volume limit. The leading piece is contained in the Bethe Ansatz 
(identical to the asymptotic Bethe ansatz), while the leading virtual correc- 
tions computed from generalized multiparticle Liischer formulas provide a 
way to compute exactly the 4-loop wrapping corrections. Thus we may split 
the 4-loop coefficient of ([1]) into 

78 (M) = 78 Bethe (M) + ^ rappin9 (M) (2) 

The Bethe ansatz term has been computed in [35J and can be found in 
table 1 of that reference. The aim of this paper is to compute the wrapping 
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part <y™ rappm9 (M) from multiparticle Liischer formulas for the AdS§ x S 5 
worldsheet QFT. This is a generalization of the computation of the 4-loop 
anomalous dimension of the Konishi operator [28] which corresponds to 

7 ™ m ™(2) = 324 + 864C(3) - 1440C(5) (3) 

The function A(M, J) has numerous interesting properties. Firstly, its 
large M limit is related to the cusp anomalous dimension [48] : 

lim A(M, J) - J - M ~ 2 lcusp (g) logM (4) 

M— >CJO 

It is therefore expected that the coefficient of log M does not depend on J. 
The cusp anomalous dimension can be investigated both from the perturba- 
tive side [49] and from the strong coupling side [501 ED EH E2J [53j EU EH] with 
an interpolating answer following from the BES equation derived from the 
asymptotic Bethe ansatz |1U| . The applicability of the Bethe ansatz analysis 
rests on the independence of '"fcusp(g) on J so that wrapping interactions do 
not contribute in this limit. This can be argued both on the perturbative side 
[56] and on the strong coupling string side [57], but it would be interesting 
to verify directly that wrapping contributions will vanish in this limit. This 
is among the aims of the present paper where we verify this for J = 2 and 
4-loop order. Once these wrapping corrections are obtained, thus completing 
our knowledge of 4 loop anomalous dimensions of twist two operators, there 
are also other very interesting features of their large M limit which could be 
investigated in particular its 'reciprocity' properties [58j, further studied in 

Secondly, the coefficients of (PQ) at /-loop, 72/ (M) are expected to obey, 
for M = 4 SYM, the principle of maximum transcendentality [21]. This 
means that they are expressed in terms of (nested) harmonic sums and ( 
functions such that the degree of transcendentality of 72/ (M) equals 2M — 1. 
The degree of transcendentality for a product is defined to be the sum of 
the degrees of each factor. The degree of transcendentality of ((n) is defined 
to be n while the degree of the harmonic sum 5'j lj j 2v .. i j fe (M) is £V \ik\- The 
harmonic sum «S'i 1 ,t 2) ...,i fc (M) is defined recursively as 

5* llfal ... A (M)=X; J ^ig^5ia,..A(n) (5) 

n=l 
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while the elementary harmonic sums with a single index are given by 

sign (.7) 



n=l 

In [35] it was verified that the part of the 4-loop result following from the 
asymptotic Bethe ansatz indeed is composed of terms of transcendentality 
degree 7. We would like to verify in the present paper that the contribution 
of wrapping corrections will also have the same degree of transcendentality. 

The third and perhaps the most nontrivial property of the anomalous 
dimensions (CQ) is the interrelation with BFKL |36j and NLO BFKL [60] 
equations which describe perturbatively gauge theory high energy scattering 
in the Regge limit. From this formalism one obtains a specific prescription 
for the pole structure of the analytical continuation of A(M, J = 2) to M = 
— 1+uj. The /-loop coefficients of (pQ) have to have the following pole structure 
around M = -1 + u [5BJ: 

T2 (u,) ~ -4(2+0 + O(a;)) (7) 



7 4 (w) 



. _ 16 (£ + £ + oK) ) (8) 

~ -64^ + ^ + (I)) (9) 

~ ^(^ + M + (y) (1 0, 

where the leading pole comes from the BFKL equation while the second one 
is a consequence of NLO BFKL. It has been verified [35] that the l-,2- and 3- 
loop result exactly agrees with the above BFKL and NLO BFKL predictions. 
The main conclusion of [35] was that the Bethe ansatz part of 7s does not 
satisfy this constraint. Its expansion around M = — 1 + uj is3 

.tf-V) „ 256 ( 4 + A + HP) _ _ + f J.^ 

1 cj 7 a; 4 u 3 \uj 2 1 1 



The main motivation of this work is to compute the corresponding wrapping 
contribution ^ rappmff (M) and to check whether the sum 7s(co>) = 7^ e * /ie (co | ) + 



1 We are grateful to Vitaly Velizhanin and the authors of 35] for informing us of this 
explicit expansion. 
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^wrappmg^ ^ correc ^- ana lytical properties required by LO and NLO 
BFKL. 
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3 Wrapping correction for twist two opera- 
tors 




Figure 1: Multiparticle Liischer correction. The vertical lines represent the 
physical particles forming the multiparticle state, while the double line loop 
represents the on-shell 'virtual' particle with complex momentum. 

The aim of this section is to extend the calculation of the leading wrap- 
ping correction of the Konishi operator to generic twist two operators. This 
amounts to calculating the leading finite size energy correction of a specific 
M particle state. In [28] the authors conjectured a Liischer type formula 
for generic multiparticle states and successfully applied it for the Konishi 
case: M = 2. The formula consists of two parts: The first describes the ef- 
fect how the finite volume modifies the particles' quantization conditions and 
shifts their momenta, while the second is due to virtual particles propagating 
around the cylinder and directly changes the energy as 

AE(L) = ~J2 / STr - [S a al a a^Pi)S a a ? a {q,P2) • • • S^Mpm)] e~ ? - (9)i 

(12) 
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The content of this formula is schematically represented in Fig. 1. This 
formula applies for an M particle state with all particles of type a, such 
that their consecutive self-scatterings preserve this state and determine their 
momenta pi by the BA equations. The matrix S^(q,p) describes how a 
virtual particle of type b with momentum q scatters on the real particle of 
type a and momentum p. The exponential factor can be interpreted as the 
propagator of the virtual particle. What makes difficult to apply this formula 
in practice is that we have to sum over all possible virtual particles in the 
theory (both fundamental magnons Q = 1 and the infinite tower of their 
bound states Q > 1), their polarizations a\ and, even more, over all possible 
intermediate states a 2 , . . . , a at- 

Let us focus on the leading wrapping correction for the anomalous di- 
mension of a twist two operator. Since the modification of the BA equation 
appears at order g 8 its contribution to the energy will be subleading and 
it is sufficient to analyze equation ffl2l) . In the next section we explain all 
the ingredients of this equation specified to the twist two case, while in the 
subsequent one we perform the actual calculation. 

4 Main ingredients of the wrapping correc- 
tion 

In order to apply formula ffl2|) for twist two operators we have to explain 
the following ingredients: What is the state such an operator corresponds 
to and how are their momenta determined by the BA equation? What are 
the virtual excitations, what is their exponential damping factor and how 
they scatter on the state that corresponds to the twist two operator? Let us 
investigate these questions in this order. 

Multiparticle states corresponding to twist two opera- 
tors 

The state which corresponds to twist two operators is an M particle state 
with M even. According to the sw(2|2) (g) sii(2|2) classification each particle 
belong to the si (2) sector of the fundamental representation, which is re- 
alized in terms of the totally anti-symmetric representation, thus has label 
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a = (1, l)o The volume is L = 2 and the momenta of the particles are 
determined by the BA equation. At leading order in g their rapidities 



1 p I j) 
u(p) = - cot - * / 1 + 16g 2 sin 2 - 



are given as the roots of Baxter's Q function, Pm{u). This is a polynomial 
of order M which is given explicitly (at 1-loop levejf]) by the generalized 
hypergeometric function (see [62], [63] and the methods of [M] ) 

Pm(u) = 3 F 2 (-M, M + 1, i - iu; 1, 1|1) 
Since this polynomial is even each rapidity Uk comes in pairs: 

M 

Pm{u) oc — Uk) = (u — Ui){u + Ui) . . . (u — Um)(u + Um) (13) 

k=l 

This is a useful form that we will use later on. 



Virtual particles and the exponential factor 

The virtual particles turn out to belong to the completely anti-symmetric 
representation of su(2\2) ®su{2\2) first discussed in [65]. This representation 
exists for any integer Q and has dimension (4Q) 2 . The dispersion relation of 
their particles leads at leading order to the exponential factor 

e -e Qiq) L = e -2La rcsi nh^ = ( Z _Y = ^ ^ 

\Z+) "(g2 + Q2)L 

where here and later on we use that 

Observe that only z + scales like g~ l but z~ goes as g in the weakly coupled 
regime. 

2 We use the convention in which fermionic coordinates have labels 1,2 while bosonic 
ones 3,4. See [28] for details. 

3 This expression has been generalized to higher loops in [61] . However we will not need 
these expressions for our computation as they would contribute only at higher loop orders. 
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The scattering matrix 



In [28] it was described how a particle with charge Q and parameters z ± 
scatters on a fundamental particle with Q — 1 and parameters x ± 

±, ^ 2u±i (, L 16 9 2 \ 

X ± (U) = 1 + A 1 - T7T 

K } 4<? \ V ( 2u±i ) 2 J 

The scattering matrix can be determined, using the superfield methods of 
[66], from the requirement that it commutes with the symmetry charges of 
su(2\2) ® su(2\2) up to a scalar factor which was obtained from the bootstrap. 
Let us recall the scalar and matrix part of the scattering matrix. 

Scalar part 

In the bootstrap procedure the composite particle z is realized in terms of 
its individual constituents zf such that the bound-state condition is realized: 
(z~ = Zi,zf), {z^ = zf,z+ = z^),...,{zq = z+_ llZ + = z+). The scalar 
factor of the Q— 1 scattering is then nothing but the product of the individual 
scalar factors of the 1 — 1 scatterings: 

8=1 

where 

s?^(z±,x±) = z ~ + - x+ _ 1 -^r (15) 

Z + — X 1 

The choice for the constituents is not unique and such configuration was 
adopted which had the most zf parameters of order g^ 1 . This leads to the 
result which, in terms of q and u, can be written as: 



*J(2), A (q-i(Q-l)-2u)(2u + iy 



Matrix part 



(q - i(Q + 1) - 2u)(q + i(Q - 1) - 2u)(q + i(Q + 1) - 2u) 



In describing the matrix part we can exploit the fact that it is a tensor product 
of two copies of the same S-matrices. This reflects the su(2\2)®su(2\2) nature 
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SB(q,u,j + l) 



of the symmetry. Thus it is enough to analyze one copy of the anti-symmetric 
representations. It has 4Q states: 2Q bosons and 2Q fermions. From the 
sm(2|2) symmetry it is possible to calculate how these particles scatter on 
the particles of the fundamental representation, see (28] for the result. 

For calculating the wrapping correction we need only the Sjj , I, J — 
1, . . .AQ part of the matrix. We will interpret it as a matrix acting on the 
anti-symmetric representation of charge Q and denote it by Sq!^ atrix (u, q) . 
Lets us go through systematically all possible matrix elements. We use the 
notation of [28J adapted to anti-symmetric representations. The first bosonic 
state 1 = 1 can scatter into itself 

SB0(q,u) :=SH(q,u) 

The bosonic state I = j, j = 2, . . . ,Q can scatter either into itself, or into 
the other bosonic state J = j + Q. In a similar manner this other bosonic 
state with J = j + Q can either scatter into itself or scatter back to / = j. 
Thus for each j = 2, . . . , Q these scatterings form a 2 x 2 matrix 

SH(q,u,j) Slj +Q (q,u,j) 

The boson with index I = Q + 1 scatters again into itself 

SBQ(q,u) = Sl^l(q,u) 

The fermionic part of the scattering is much simpler, since it is diagonal and 
is the same for both fermionic particles: 

SF(q,u,j + 1) = Sjgg(g,«) = SH%+i(q,u) 

The explicit form of these matrix elements can be found in the Appendix of 
|28j . Here, for convenience, we grouped them slightly differently by shifting 
some of the labels. (For the paper to be self contained, we list the needed 
matrix elements explicitly in Appendix D). 

5 The calculation of the wrapping correction 

Having introduced all the needed quantities we can formulate the wrapping 
correction as 

M M 



* e = -e f I n^.^^n 3£«.<«. "-)) v ! ""'' 

q = 1 J -00 k = l k = 1 



QmatrixVii k) ) J / 2 _|_ Q 2 y 
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Where the super-trace part can be further elaborated as 



M MM 
? *J(2) 



,u k ) + ,u k ) 

k=l k=l k=l 

Q-1 M Q-1 / M \ 

2 e n sf ^ «*> + e Tr n ^> ^> ( ie ) 

i=0 fc=l j=l \fc=l / 

We evaluate this part carefully in Appendix A and here we just summarize the 
main steps. Since we expect the wrapping contribution to appear at order 
g 8 , and there is an explicit g 4 dependence coming from the 'exponential' 
part, the super-trace part has to be of order g 2 , so its leading term should 
vanish. We check this requirement first and then make a systematical small 
g 2 expansion of each term later. Schematically we expand each function 
as /(<?) — /o(l + 9 2 3f + 0(g 4 )). The key feature why we were able to 
do the calculation is, that the leading order part of the bosonic matrices 
SB(q,u,k ,j)o can be diagonalized in a w^-independent way. This enables us 
to rewrite the bosonic matrix contribution as a sum over the contributions 
of the two eigenvalues SB ± (q,Uk, j)o- Explicitly we found that 

. , N i + 2ij + q — iQ — 2u 



SB + (q,u k J)o 



i + q — iQ — 2u 



Additionally, when we extended the summation for these eigenvalues (from 1 
to for + and form Q — 1 to Q for — ) we could incorporate the contributions 
of the two separate bosonic terms SBO(q,Uk)o and SBQ(q,Uk)o- We also 
observed that SB~ (q,Uk, j + 1) = SB + (q, w fc , jOf^zf such that the — bosonic 
summation can be also shifted and the zero order part for wrapping has the 
form: 

m Q-i M Q-1 M 

(n SSL,*?, «o) = -2 e n sf ^ u ^ e n sb+ ^ **> ^ 

k=l j=0 fc=l j=0 k=l 

Further checking that the fermionic parts and the bosonic parts are the same 
for each j, (SF(q,u k ,j) = SB + (q,u k J)^^r), we could see that they 
completely annihilate each other at this zeroth order. 

In doing the calculation at first order we performed the same steps: We 
diagonalized the bosonic contributions upto 0(g 2 ), we extended the bosonic 
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summation to incorporate the extra separate bosonic pieces, shifted the sum- 
mation for the second bosonic part, and finally exploited the fact that the 
fermionic and bosonic zero order terms are the same for any j. As a result 
we arrived at formula 

M Q-l M M 

k=l j=0 k=l k=l 

where after some lengthy but straightforward calculations we found that 

1 1 
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5SBF(q,u k ,j) = 2 
1 + 4«r 



2j-iq-Q 2(j + l)-iq-Q 



Exploiting the very simple u k dependence of the summand we can recognize 
the one loop BA energy and replace it with the harmonic sum S\(M) as 



M M 

In order to abbreviate future formulas we will suppress the arguments of 
harmonic sums if they are M. The equation (fT3i) can be used to write the 
result in a more economical way: 

M SB (q, Uk,J)o = — — -jY- — -tt^— 

H PM(m - i(Q - 1))) 

The final form of the wrapping correction 

We arrived at one of the main results of the paper. We merely have to 
collect all the ingredients of the wrapping correction: the exponential factor, 
the scalar factor and the matrix part. There is an elegant way of writing the 
whole wrapping correction in terms of Baxter's Q function^ as 

.„ „ lc! f r dg T M (q,Q) 2 16 

= -«« «. t, r mM) wrcrr (17) 



4 We denote the Baxter's function as Pm(-) in order to avoid confusion with the Q 
parameter. 
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where 



i?M(g,g) = PM(i(g-^(g-i)))PM(^(g+^g-i)))PAf(i(g+«(g+i)))PM(i(g-*(g+i))) 



and 



Q-i 



T M (q,Q) = Yl 

j=0 



PM\^{q-i{Q-l))+ij 



2j -iq-Q 2(j + l)-iq-Q\ \ 2 
An alternative way of writing this function is 

TmM) = ^(^-w-i))) _ ^(^ + -w-D)) +f . ffa|Q) 



q-iQ 



q + iQ 



where 



J Pm (|(<Z - i(Q - 1)) + ij) - Pm (\(q - i(Q + 1)) + ij) 



It is not difficult to show that TM(q,Q) defines a polynomial. 



Wrapping correction for odd particle number 

In computing the anomalous dimension of twist two operators the consider- 
ation of odd states played an important role [35J. This meant an analytical 
continuation of the result from even to odd M. The hypergeometric function 
determines the allowed rapidities: they are all paired except one which has 
u = 0. This is not physical since it corresponds to p = it and not p = 0. 
Nevertheless, we can try to extend the wrapping correction to this case. The 
first problem immediately arises at zeroth order. The extra p = n state 
gives different extra factors to the two bosonic and fermionic contributions. 
Namely the B + part is unchanged, the B~ gets an additional —1 and annihi- 
lates B + , while both fermions got an extra i. Since we would like to maintain 
for the analytical continuation that the correction starts at g 8 we suppose 
that in a theory where u = is an allowed state the symmetry between the 
two fermions are broken and their contributions have opposite signs. Thus 
we add extra factors as i and —i to the two fermionic components. With this 
assumption we could calculate the wrapping correction of those odd states. 
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The only difference compared to the even one resides in the function T(q, Q) 
which turns out to be 



T° M dd (q,Q) = J2 

3=0 



1 



1 



[2j-iq-Q 2(j + l)-iq-Q\ 



So in essence a 'nice' analytical continuation of (1171) amounts to inserting an 
additional (— 1) M sign in the definition of T M (q,Q). 



Calculating the integral 

We calculate the integral by closing the contour on the upper half plane and 
taking residues over the poles of the integrand. We have two types of poles. 
The pole at q = iQ is of kinematical origin since it does not depend on the 
scattering matrix. In contrast, the poles of the function R(q, Q) depends 
on Uk thus are determined by the dynamics and come from the scattering 
matrix. They correspond to \i terms in the Liischer correction. In the paper 
[28] an argument, based on kinematical considerations, was given that such 
terms are absent in the weakly coupled (g — > 0) limit, so we expect their 
contribution to vanish. Although they do not vanish separately but when 
we sum up their contributions for all the bound-states (for Q) the result 
is indeed zero. We checked this explicitly for the first few M cases. As a 
consequence in calculating the integral for the wrapping correction we can 
take the residue only at the kinematical pole q = iQ. Subsequently we have 
to sum the resulting expression for Q running from 1 to oo. 



6 Determination of the final form of ^ rappm9 (M) 

Assuming the maximum transcendentality principle, we expect that ^ rappm9 (M) 
has the following structure: 

7 ™ rappm9 (M) = C 7 {M) + C 4 (M)C(3) + C 2 (M)C(5) (18) 

where the coefficients C n (M) have transcendentality degree n. 

In Appendix B we have analytically derived the coefficients of £(3) and 
C(5) and found that indeed they have the expected transcendentality struc- 
tur<H 

C 2 (M) = C 4 (M) = -512£%S-2 (19) 

5 In the following all harmonic sums are evaluated at M so we suppress the argument. 
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It remains to determine the rational part Cj(M). A-priori the number of 
independent harmonic sums of transcendentality 7 is quite large, but we 
may use the structure of the Liischer correction i.e. the fact that Sf gets 
automatically factored out to significantly simplify the analysis. Hence we 
may write Cj{M) = SfCa(M). We are thus left with harmonic sums of 
transcendentality 5. Assuming that the index —1 does not appear we are left 
with a basis of 41 harmonic sums. 

We then computed analytically the residues of the integrand at q = iQ 
and summed up the resulting expression from Q — 1 to oo. We did this 
for M — 1 to M — 41, where we included, similarly as in [35], also the 
unphysical odd values of M. A justification of the precise form of analytical 
continuation of the even M integrands to odd M is the agreement of the 
coefficients of ((3) and £(5) and the fact that the remaining term was just a 
rational number. The results from M = 1 to M = 41 fixed all the coefficients 
of the 41 harmonic sums which turned out to be simple integers. The result 
for the rational part obtained in this way is 

C 7 (M) = 2566\ 2 {-S 5 + SL 5 + 25-4,1 - 2S 3) _ 2 + 2S„ 2 ,_ 3 - 4SL 2 ,_ 2)1 ) (20) 

As a check we then computed the wrapping correction from Liischer formulas 
for M = 42 and M = 44 and got perfect agreement with the above formulas 

(HHD-flUD- 

In the remaining part of the paper we will first verify that ^ rappin9 (M) 
does not have any contribution of the order of log M in the large M limit 
so that the cusp anomalous dimension remains unchanged. Then we will 
analyze the analytical continuation of >y™ rappm9 (M) to M = — 1 + u and 
check compatibility with LO and NLO BFKL. 

7 Large M asymptotics of <y™ rapping (M) 

Wrapping corrections should not change the leading asymptotic behaviour of 
the anomalous dimension of twist-2 operators. On the other hand the final 
formula which we obtained is of the form 

^ rappin9 (M) = 128Sf{M) x (finite when M — oo) 

where Sf(M) scales as log 2 (M) when M — >• oo. The only option for the 
wrapping correction not to change the leading large M asymptotics is that 
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the finite part vanishes at infinity. If this is the case then the first wrapping 

M 



contribution to the asymptotic behaviour may enter at the order log which 



is sub leading comparing with log(M). 

Let us check if the leading expansion around M = +00 of our result 
vanishes. The values at infinity of the nested harmonic sums can be expressed 
in terms of multivariate zeta functions (Zagier-Euler sums) as was shown in 
|67j . These sums can be reexpressed in terms of ordinary Euler sums. The 
relations between the former and the later can be found using EZ-Face — an 
on-line calculator for Euler sums [68]. The results relevant for us are 

S-a(«0 = 4 C(2) 

£-6(00) = -H C ( 5 ) 

S 5 (oo) = C(5) 
S M (oo) = -C(2)C(3) + 3C(5) 

SsMoo) = i C(2)C(3) - ^C(5) 
^-2,-3(00) = yC(2)C(3) - ^C(5) 
5-2,-2,i(oo) = I |C(2)C(3) - ||C(5) 
Plugging them into the coefficient of Sf in the wrapping correction 
-5C(5) - 4S1 2 C(3) - 2S 5 + 2S , _ 5 + 4S , 4 ,i - 4S 3 ,-2 + 4S , _ 2 ,_ 3 - 8S_ 2 ,_ 2jl 

we obtain that it vanishes as we expected. It would be very interesting to 
analyze the large M limit in more detail along the lines of [581 ES] . 



8 Analytical continuation of ^ rappin9 [M) 

The way how harmonic sums can be analytically continued was nicely ex- 
plained in [67]. There is a pronounced difference between the cases when 
all indexes are positive and when we have at least one negative index. The 
former case can be described systematically, while the latter one requires a 
case by case study. We relegate the details of the analytical continuation of 
the needed harmonic sums to Appendix C. Using the results there we can 
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list the analytical continuation of all harmonic sums of interest. Since we are 
concerned with the singularity around —1 only up to the third order pole, it 
is enough to expand Si to third order and keep only the singular part of all 
other harmonic sums: 



__ + wC(2 )-„ 2 C (3) + a/C(4) + ... 
1 

-i(C(2)^-C(3)^ 2 + C(4)^+...) + --- 



S x ( 


-1+u 


S 5 ( 


-1+U 


5*4,11 


-1 + u; 


5-5 ( 


-1 + u 


5*3,-2 1 


-1 + u 


5*-2,-3 l 


-1 + u 


5*— 2,-2,1 1 


-1 + u 



u 

1 



:(2C(-2) - 2((-3)u + 3((-4)u 2 + ...) + 



- — (3C(-4V - 2C(-3) + 



+ . 



or 
1 



-(-2S_ 2 ,i(oo) 



+a;(25'_3 ) i(oo) + 5_ 2 , 2 (oc) - C(-2)C(2))) + . . . 



In the above expressions ( of a negative argument denotes an alternating 
version of the ( function: 



C(-n) = E 



k=l 



(-1)' 

k n 



(21) 



which is linked to the ordinary ( function through the well-known relation 

C(-n) = (2 1 -™ - l)C(n) 

The harmonic sums at infinity can be expressed in terms of Euler-Zagier 
sums. We found the following relations useful: S*_2,i(oo) = 3) and 
2S*_3 i i(oo) + 5*_ 2 ,2(°o) = — 16 C (4)- Some of them can be proven, but some we 
obtained using the program EZ-Face [55] . To simplify the final form we also 
used that C(2) 2 = |C(4)- If we plug all these expressions into the wrapping 
correction 

^2 



128S? [-5C(5) - 45_ 2 C(3) - 2S 5 + 2S. 5 + AS 4 ,i - 45* 



3,-2 



45_ 



2,-3 



8S_ 



2,-2,1 
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we obtain the leading singularities around M = 



— 1 + to as 



7r -» ~ 256 ( A - Hffi + m + + o (L\ ) (22) 

\a; 7 ur u A Aui 6 \uj 1 ) J 

which, when combined with the Bethe Ansatz result f fTTl) 

n) „ 256 ( =1 + L t + !M _ _ M!) + ( L\ ) (23) 

agrees with LO and NLO BFKL expectations 

Finally we note that we can continue around any other negative integers 
to compare with other predictions. In particular, it is known that there is an 
uj^ 1 pole at even negative integers [35], whose coefficient is fully reproduced 
just by the Bethe Ansatz answer. We will thus check that our wrapping 
correction does not give any contribution to this pole which we analyze at 
M = — 2 + uj. The terms which can contribute at this order are S$ and S-5. 
Since the analytical continuations are S±s(M = —2 + uj) = — ^ there is no 
contribution from wrapping at this order which is consistent with the BA 
results. 

9 Conclusions 

In this paper, using integrability properties of the light-cone quantized world- 
sheet QFT of the string in AdS^, x 5* 5 , we have obtained the four- loop wrap- 
ping correction to the anomalous dimension of twist- 2 operators with spin 
M. As a first check we determined the large M asymptotic behaviour of our 
expression and concluded that it does not modify the Bethe Ansatz result 
for the cusp anomalous dimension. In contrast, the wrapping correction is 
essential for the correct behaviour of the twist-2 anomalous dimensions un- 
der analytic continuation toM = -l + w. Indeed it exactly cancels all the 
higher order poles in the Bethe Ansatz result and the remaining leading poles 
exactly agree with LO and NLO BFKL predictions. 

Let us note that despite the apparent complexity of the computation, the 
string theory calculation performed in the present paper is still much sim- 
pler than any corresponding gauge theory perturbative computation. This 
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suggests that one can use string theory methods of the AdS/CFT correspon- 
dence as an efficient calculational tool even in the weak coupling perturbative 
regime. 

In fact, it would be interesting to analyze in detail the interrelations 
between the perturbative and string theory computations as it might give a 
hint on streamlining higher loop perturbative methods. 

Moreover it would be very interesting to use the complete formulas for the 
4-loop anomalous dimensions to understand better the underlying theoretical 
structure both from the point of view of their asymptotic properties along 
the lines of [58l [59] and of the analycity structure for higher negative integers 
as discussed in [32] and further aspects of the link with BFKL like in [70] . 

Apart from these physics issues, there still remain several aspects of our 
derivation which may be improved. We have managed to do the exact cal- 
culation of the result for the terms proportional to the £(3) and £(5) (see 
Appendix B) but the derivation of the rational part was essentially based on 
the maximum transcendentality conjecture. It would be very instructive to 
obtain this result from first principles. 

Finally let us note that the exact form of the Luscher corrections may be 
a precision test that has to be satisfied by any yet-to-be constructed exact 
spectral equation (or more probably a set of coupled nonlinear integral equa- 
tions). Our computation in this paper and the very nontrivial consistency 
requirement with the analytical structure predicted from BFKL can be un- 
derstood as a check of the form of the Luscher corrections for the worldsheet 
QFT of the AdS^ x S 5 superstring. 
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A Evaluation of the integrand 



In this appendix we elaborate on the super-trace part of the integrand. Since 
we expect the wrapping contribution to appear at order g 8 the super-trace 
part has to vanish at leading order. We check this requirement in the first 
subsection and then make a systematical small g 2 expansion of each term in 
the second one. 



Zeroth order contribution 

Since the rapidities appear always in pairs (u, —u) we found it useful to com- 
bine their contributions. Using the explicit form of these functions from [28] 
(see Appendix D below) we evaluate the expressions entering the evaluation 
of the supertrace 

SB0(q,u)SB0(q, -u) = 1 
for any g. To abbreviate the formulas we introduce the following notation 

f(g)=f (l+g 2 5f + O(g 4 )) 

In this notation the single diagonal bosonic part is given by 

Qnn( ^ Qnn( ^ (q + t(Q - l)) 2 - Au 2 
SBQ(q,u) SBQ( qi -u) = {q _ l{Q _ 1))2 _ Au2 

While the diagonal fermionic part for j — 0, . . . , Q — 1 reads as 
SF(q,u,j) SF(q, -u,j) = {q _ l{Q _ l)? _ Au2 

The main complication is in the matrix part. One can observe, however, 
that evaluating the matrix part at g = the resulting matrix SB(q,u, j)o 
has eigenvectors (^few, 1) and ( q 2+Q2 ? 1) ? which are independent of the u^- 
Thus we can diagonalize all of them simultaneously leading to 



GSBiq^j^G- 1 



i+2ij+q—iQ~2u g 
i+q— iQ— 2u 

g 2ij+q-i(l+Q)-2u 2u-i 



i+q—iQ—2u 2u+i 

Let us denote the corresponding eigenvalues by SB ± (q,u, j)o- As a conse- 
quence 

(M \ M AI 

Y[SB(q,u k ,j) ) = Y[SB + (q,u k ,j) + JJ SB~ (q, u k , j) 
k=l J o k=l k=l 
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where we have 

SB (q,u,j) SB (g, = {q _ i{Q _ w _ Au2 

Now couple of observations are in order. We can see that 

SB + (q, u, 0) SB + (q, -u, 0) = SB0(q, u)SB0(q, -u) = 1 

and that 

SBQ(q, u) SBQ(q, -u) = SB~(q, u, Q) SB~(q, -u, Q) 

This means that we can incorporate the contribution of SBO into SB + by 
extending the summation over j in (1161) from j — 1 to j — 0. In a similar 
manner, the summation of the contributions of SB~ when extended naively 
to j = Q turns out to automatically incorporate the contribution of SBQ. 
We can further realize that 

SB + (q,u,j) SB + (q, -u,j) = SB~(q,u,j + l) Q SB~(q, -u,j + 1) 

Thus the two bosonic eigenvalues contribute in the same way This means 
that at leading order the super-trace part reads as 

M Q-l M Q-l M 

(n s^Satrd^ «*)) = - 2 e n sf ^ u ^ e n sb+ ^ «*. ^ 

k=l j=0 k=l j=0 k=l 

But finally we can observe that the bosonic and the fermionic contributions 
are the same for each j: 

SB + (q,uJ) SB + (q, -u,j) = SF(q,u, j) Q SF(q, -u,j) 

so the super-trace part vanishes at leading order and consequently the wrap- 
ping contribution also vanishes at order g A . 

First order contribution 

We expand now each function to the order g 2 as follows: The expansion of 
SBO is trivial, the correction vanishes. The other diagonal bosonic contribu- 
tion reads as 



M M 

]JSBQ(q,u k ) = (j[SBQ(q,u k ) 

k=l k=l 



M 

l+g 2 Y,SSBQ(q,u k ) 



k=l 
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The fermionic is given by 



M M 

]JSF(q,u k ,j) = (j[SF(q, Uk ,j) ) 



k=i 



k=i 



M 



l + g 2 Y,SSF(q,u k ,j) 



k=l 



In the matrix case we use the same matrices G we used to diagonalize 
SB(q,u, j)o and bring G SB(q,u, into the form 

SB+(q,u,j) (l + g 2 5SB+(q,u,j)) 0(g 2 ) 

0(g 2 ) SB-(q,u,j) (l+g 2 5SB~(q,u,j)) 



The advantage of this form is that 



vfc=l 



Tr[l[SB(q,u k ,j)} = (j[SB + (q,u k ,j) ^ 

M 

(j[SB-(q,u k ,j) ) 



M 



k=l 



l+g 2 J2^B + (q,u k ,j) 
k=i 

M 

l+g 2 J2^B-(q,u k ,j) 



+ 



k=i 



We have checked again by explicit calculation that 

6SB + (q, u, 0) + 5SB + (q, -u, 0) = 

and that 

5SBQ(q, u) + 5SBQ(q, -u) = 6SB~(q, u, Q) + 5SB~(q, -u, Q) 

Thus both bosonic summations can be extended as before. Shifting the 
summation in the SB~ case we can bring the whole leading correction into 
the form: 



k=i 



M Q-l M M 

j=0 k=l 



k=l 



where 



5SBF(q,u k ,j) = 5SB + (q,u k ,j)+5SB-(q,u k ,j + 1) - 25SF(q,u k ,j) 
We could explicitly calculate this quantity, which turned out to be 



5SBF(q,u k ,j) 



16 



1 + 4x4 



1 



2j-iq-Q 2(j + l)-iq-Q 
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B Exact calculation of the coefficient of £(3) 

and ("(5) 

Having taken the residue at q = iQ we will encounter derivatives of the 
Baxter's Q functions (denoted here by Pm(-) to avoid confusion with Q) 
around four different points, namely around ±| and ±^+iQ. In the following 
calculation it is more convenient to use the function with a shifted argument 

U{±Q)=P M {± l - + iQ) 

It has the following expansion around Q = : 

U(Q) = 1 + 2Sr Q + 2(S 2 + S_ 2 ) Q 2 + . . . 

After taking the residue of the integrand we have a sum of rational functions 
of Q where the denominators are products of powers of Q and U(Q)U(—Q). 
We can make a partial fraction expansion of this result and separate the 
terms having denominator Q n : 

^ -44 -43 -42 ^1 ^ 5 AjjjQ) 
Q 5 Q A Q 3 Q 2 Q ^ U(QYU(-Qy 

Explicit calculation showed that the quantities of interest read as 

A 5 = 10 ; A 3 = U"(0) - 4U'(0) 2 + 12(d Q T(q, Q))\ q=lQ=0 

In writing A% into the form above we used U(0) = 1 and the previously cal- 
culated values of A 5 and At = 12(T(0, 0) + U'(0)) = 0. The other coefficients 
A 2 and A x and especially the polynomials Aij(Q) are very cumbersome to 
compute. Using the expansion of U(Q) in terms of harmonic sums together 
with {pQT(q, Q))\ q=i Q =Q — Sf — S , _ 2 we obtain the form of the wrapping 
correction 

AE = -128g 8 Si [5C(5) + 45_ 2 C(3) + . . . ] 

where the ellipsis denotes a rational number which can be expanded in terms 
of harmonic sums and which we determine in section 6. 

We note that using perturbative methods in the M = 4 super Yang-Mills 
theory the coefficient of £(5) was calculated in [69]. It is exactly reproduced 
by our computation. The conjecture for the coefficient of £(3) given in [69] 
disagrees, however, with our result. 
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C Analytical continuation of harmonic sums 



There is a pronounced difference in the analytical continuation of harmonic 
sums between the cases when all indexes are positive and when we have at 
least one negative index. The former case can be described systematically, 
while the latter one requires a case by case study. 



Analytical continuation with only positive indices 

The analytical continuation of harmonic sums with all indices positive can 
be done inductively. One starts with the simplest one S a (n) and uses the 
general strategy to move the variable n from the upper bound of the sum to 
the summand as 

(OO OO \ -. OO 

E-E r *)-E^ 
j=l j=n+l/ J k=l K > 

Since we are interested in the analytical continuation around — 1 we explicitly 
separate the singular and regular pieces as 

S a (-1 + x) = ~ + S a {oo) - d a (x) 
x a 

where we found it useful to introduce the function 

This function is regular around x = and has the expansion 

d a (x) = ((a)-xa((a+l)+x 2 + ^(a+2)+. • -+(-l) n x n (° + ^ ~ ^ C(a+n) 

Suppose now that we have already analytically continued ^...^(n) and then 
we want to continue analitically S a ^,^ c {n). Using the previous strategy we 
can write 

(\ | x i 

J2~ - a s K-,c{3) = S aA ...Aoo)-J2 (k T )a Sb,...,c{k+n) 

j=l j=n+l/ J k=l ^ ' 
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Let us specify this formula for the case of interest, namely to S at b(n) around 
n = — 1. Focusing only on the singular part, which comes from the k = 1 
term, we have 

S a ,b(~ l +x) = — - [S b (x)] + reg 
x a 

where by "reg" we mean terms regular for x — > and 

S b {x) = S b {oo) - d b {x) = xb({b + 1) - x 2 f 6 + 2 *J C(6 + 2) + • • • 
This is all we need for the analytical continuation for positive indexes. 

Analytical continuation with at least one negative index 

In the case when we have at least one negative index in the harmonic sum we 
have to be careful. The main problem is that the harmonic sum determines 
two different analytical functions, depending on whether we continue from 
even or from odd values. Since in our problem we continue from even ones 
we have to extend our original harmonic sum to odd values first, then to do 
the analytical continuation. Since the way how one can do this extension 
for the functions we need is thoroughly explained in |67J we merely cite the 
result. For instance, in the simplest case when we have just one index the 
function 

S+ a (n) = (-l) n £- (n) + (1 - (-l)»)£_ a (oo) 

is the same as S^ a (n) for even values and is extended to odd values in such 
a way that it can be described by one function. In calculating the analytical 
continuation around —1 we use the method of moving n from upper bound 
to the summand and obtain 

St a (-1 + x) = ^ + S_ a (oc) + d_ a (x) 
x a 

where the function d- a (x) is regular around zero and have the expansion: 

cL a (x) = C(-a) - xaC(-a - 1) + x 2 f ° ^ J ((-a - 2) + . . . 

Similarly one can define the extension of the harmonic sums with the first 
index being negative to be 

= (-l) n ^,...,W + (1 - (-l)")SL a , 6 ,... c (00) 
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as the proper continuation from even values to odd ones. Its analytical 
continuation around —1 has the singular part 

S\ b ,..J-l+x) = ^- a [S b ,..., c (x)]+Teg 

The analytical continuation of the function S a) -b{n) is more tricky and one 
has to use the functional relations the harmonic sums satisfy to define 

= (-l) n 5 a ,_ fe (n) + (1 - (-l) n )(S a ,_ 6 (oo) - S b (oo)(S a (oo) - S a {n)) 

Thus the singular part around —1 reads as 

We also need the analytical continuation of S-2,-2,i( n ) so we need 

n 

St a ,- b , c (n) = E ~S\cU) + 5_ 6)C (oo)(4+(n) - S a (n)) 

3=1 3 

Its analytical continuation has a singular part: 

^a,- 6) o(-l +*) = —S^W + ^,c(00)(^ - ^) 

in which we have 

S\ c (x) = S^ c (oo)-J2jj^ b Sc(k + x) 

= x(bS-b-i, c (oo) + c5_ bjC+ i(oo) - cS- b (oo)S c+ i(oo)) + . . . 

D Explicit S-matrix coefficients 

In this appendix we collect the explicit expressions for the matrix part of the 
Q — 1 scattering matrix derived in [28] . We choose the normalization as 

SB0(q,u) = a\(q,u) = 1 

The bosonic matrix part reads as 
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2g=*aJ(g, u) + ia 2 2 {q, u) 3 ^E^al{q, u) 

where j = 1, . . . ,Q — 1. The last bosonic part is simply 



SQ{q,u) = gq77 a l(<?,w) + ^a^u) 
while the fermionic part is given by 

SF(q,u,j) = Q—j-a 6 6 (q,u) + 3 -a 7 7 (q,u) 
The various matrix elements read in terms of x ± (u) and z (q) as 

6 r\ 

a 6 = Q— -\ — ; a 



6 —2; /x+ _ 7 2 z (a; — z + )(l — x z + ) /x H 

x + — z~ V x~ ' Q z + (x + — z~)(l — x~z~) V or 



a 10 
u io 



2 z-(x--^ + )(l-x + z + ) 
Q-lz+(x + -z-)(l-x~z-) 



2 .Q — 1 z (x — x + ) 4 .Q — 1 (z — z + ) 2 (x — x + ) 

a ± = ~ l Q Z +(X+-Z-)(1-X-Z-) ' a2=% Q (X+-Z-)(1-X-Z-) 

a 2 = ~ nn \n\ +( + ^T71 [2z~z + {Q + x~z~ - (1 + Q)ar* + ) 

Q(l + Q) — z ){l — x z ) L 

+2x + (z + + z~(-l + Q(-l + 



4 2Qx+x-(x+ - z-)(l - L vvv 7 v 



(2 — x + z )) — (x — x + )x + z [z — z + ))~j 
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